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Abstract 



.We study the Ughtest masses in the fermionic sector of an anomalous U(iy extension of the minimal supersymmetric 
Standard model inspired by brane constructions. The LSP of this model is an XWIMP (extremely weak interaction 
' } particle) which is shown to have a relic density satisfying WMAP data. This computation is carried out numerically 
after having adapted the DarkSUSY package to our case. 
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1. Introduction 

There has been much work recently to conceive an in- 
.tersecting brane model with the gauge and matter content 
of the Standard Model (SM) of particle physics [i[3,|l,ll. 
.One of the features of such models is the presence of ex- 
tra anomalous gauge I7(l)'s whose anomaly is cancelled 
via the Green-Schwarz mechanism. This characterizes this 
class of models with respect to those which still have ex- 
|tra U{1) gauge symmetries but cancel the anomaly as in 
the SM or in the minimal supersymmetric SM (MSSM) 
.(see [5i] for a review). The latter are inspired by (super) 
grand unified models or by plain extensions of the (MS)SM 
and their quantum numbers with respect to the new gauge 
simmetries are fixed by the condition of setting to zero all 
.the triangle diagrams in which at least one vertex of such 
diagrams has a C/(l)' current coming from the extra gauge 
jsimmetries. On the contrary, in brane inspired models, 
these quantum numbers are not fixed. We have introduced 
and discussed the signatures of this model in and @. 
Another important issue which deserves careful scrutiny is 
|the compatibility of this model with the WMAP data: this 
was done in [6^ (see also for related work) . The mass 

.matrix of the fermions uncharged with respect to the U(l) 
of the gauge group of the SM, gets now two new contri- 
butions in this model: one coming from the superpartner 
of the Stiickelberg boson (Stiickelino) and the other from 
the superpartner of the gauge boson mediating the extra 
[/(!)'. By taking some simplifying and reasonable assump- 
tions (to be detailed later on) on the fermion masses enter- 
ing the soft supersymmetry breaking lagrangian, the LSP 
turns out to be the Stiickelino. It is also easy to realize 
that, given the simplifying assumptions mentioned above, 
the LSP is interacting with the MSSM particles with a 
coupling suppressed by the inverse of the mass of the ex- 
tra gauge boson of the theory which must be at least of 
the order of the TeV for phenomenological reasons. The 
LSP is then an XWIMP (Extra Weakly Interacting Mas- 
sive Particle) , a class of particles which have already been 
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12|. The cross section of these 
LSPs is too weak to give the right relic abundance. This 
is why one has to resort to coannihilations with NLSPs. 
In our case, the cross section for the annihilations of the 
LSP with the NLSP and that for the coannihilations of 
the two species differ for some orders of magnitude. Once 
again this situation is not new in literature (isl 14 , 1^ but 
needs to be treated carefully: the two species will not de- 
couple as far as there will be some MSSM particles to keep 
them in equilibrium. Moreover these particles must be rel- 
ativistic so that their abundance is enough to foster the re- 
action. These points were assumed in our previous paper 
on this subject Q to allow for a simplified treatment of 
the relic abundance. In this paper we drop all simplifying 
assumptions and discuss the most general case, in which 
the extra MSSM sector is not decoupled to the MSSM 
sector and the LSP is mainly a mixture of Stiickelino and 
primcino, with small MSSM contribution. We solve nu- 
merically the Boltzmann equation in this general case. To 
do this we have modified the DarkSUSY [l3| package to 
keep in account the new interactions typical of our model. 
The results we obtain are qualitatively compatible with 
the findings in [6|: there is an ample region of the param- 
eter space which leads to a relic density compatible with 
the WMAP data. This is the plan of the paper: in section 
[2] we describe the neutral mixing in our model and sketch 
the way in which it affects the interactions. In section [3] 
we describe how we change the DarkSUSY package and 
compute the relic density. 



2. Neutral mixing 

Our model Q is an extension of the MSSM with an 
extra U{1). The charges of the matter fields with respect 
to the symmetry groups are given in table 1. 

The anomalies induced by this extension are cancelled 
by the GS mechanism. Each anomalous triangle diagram is 
parametrized by a coefficient hif'' (entering the lagrangian) 
with the assigment: 
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Table 1: Charge assignment. 
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where Ms and Mq are the soft masses of the stiickelino 
and of the primeino, respectively and u„, Vd are the vevs 
of the Higgs fields. 
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The mass of the extra boson is parametrized by Afy(o) = 
46350, where go is the coupling of the extra U{iy. The 
terms of the Lagrangian that will contribute to our calcu- 
lation are @, 0] : 



-'Stiickelino — 



'-^s^^d^i^s - V2&3^sA(") 



2V2 ^ 



^XWa'^a''FlPJi,s + (0 o 1) 



(6) 
/i.e. 



As already said, we want to deal with the case of gen- 
eral neutral mixing. The neutral mixing matrix is: 



(7) 




o 2 

Defining an = go Qh^ ^u'^ have at tree level: 



cw -sw sw\ Jg\ +gi an 

M = I Sw Cw Cw \J g\ + .92 

cwg2 + Sw gi) an 1 

(8) 

where cw is cos{9w), sw is sin{9w)- gi, 52 are the cou- 
plings of the SM electro-weak SU{2) x U{1) group. The 
structure of this matrix leaves the electromagnetic sector 
and the related quantum numbers unchanged with respect 
to the MSSM ones. 

Last, we remember 0] the general form of the ncutrali- 



3. Numerical computations and results 

Following our previous work [6] , in which we have stud- 
ied separately the case in which the NLSP is a bino-higgsino 
from that in which it is a wino-higgsino, we performed this 
general study in the same way. The plot we will show are 
generated with a modified version of the DarkSUSY pack- 
age, in which we added the new fields and interactions 
introduced by the anomalous extension. 
The free parameters that we use in our numerical simula- 
tions are the seven ones used in the MSSM-7 model: the n 
mass, the wino soft mass M2, the parity-odd Higgs mass 
Maq, tg(3, the sfermion mass scale ms/, the two Yukawas 
at and at- We add to this set five parameters which de- 
fine the U{iy extension: the stiickelino soft mass Ms, the 
primeino soft mass Mq, the C/(l)' charges Qh^, Qq, Ql- 
As our actual aim is to study the model without any sim- 
plification, we can't have control on the mass gap between 
the LSP and NLSP, because fixing it would require a con- 
straint on the masses. So we choose to let all parameters 
unconstrained and therefore to collect data in the mass gap 
ranges ^ 5%, 5 10%. In each case we have started our 
study scanning the parameter space in search of the region 
permitted by the experimental and theoretical constraints, 
i.e. the region in which we could satisfy the WMAP data 
with a certain choice of the model parameters. After that 
we have numerically explored this regions to find sets of 
parameters that satisfy the WMAP data for the relic den- 
sity. We have found many suitable combinations for both 
types of NLSP. So we have chosen some of these succesful 
models and we have computed the relic density keeping 
constant all but two parameters and plotting the results. 
We have found regions of the parameter space in which 
the WMAP data are satisfied for mass gap over 20%, but 
in the following we will only show results for the regions 
5% and 5 10%, because they are more significant. 
For simplicity we will refer to these regions as "5% region 
"and "10% region" respectively. 



3.1. General results 

In this section we want to list some results that are 
valid for both types of NLSP. 



2 



First of all^eiven the constraints on the neutral mixing de- 
scribed in [5|, we have obtained that — 1 < Qhu 1- This 
imphes that also in our general case the mixing between 
our anomalous LSP and the NLSP is small. 
We have checked that there are suitable parameter space 
regions in which the WMAP data are satisfied and we have 
found that this is true for all possible composition of our 
anomalous LSP. 

We have also checked that in each region we have studied 
there are no divergences or unstable behaviours in our nu- 
merical results. 

We have verified that the relic density is strongly depen- 
dent on the LSP and NLSP masses and composition while 
it is much less dependent on the other variables. Anyway 
it can be shown that there are cases in which the parame- 
ters not related to the LSP or NLSP can play an important 
role. We will show an example of this case in a forthcoming 
subsection. 

3.2. Bino-higgsino NLSP 

If the NLSP is mostly a bino-higgsino we have a two 
particle coannihilation. 

We chose two sample models which satisfy the WMAP 
data [iBl with mass gaps 5% and 10%. We study these 
models to show the dependence of the relic density from 
the LSP composition and from the mass gap. To obtain 
these result, we have performed a numerical simulation 
in which we vary only the stiickelino and the primeino 
soft masses. The results are showed in figure 1, with the 
conventions: 

• Inside the continous lines we have the region in which 
(^h'^)wMAP ^ ^h'^ 

• Inside the thick lines we have the region in which 
(flh'^ ~ ■ia)wMAP < ^h"^ < i^h'^ + '^o-)wMAP 

• Inside the dashed lines we have the region in which 
(n/i^ — 5a)wMAP < ^h'^ < {^h^ + 5cr)wMAP 

• Inside the dotted lines we have the region in which 
(fi/l^ - lO(7)wMAP < ^ti^ < (^^/^^ + 10o-)wMAP 

Going from the region with a 5% mass gap to that with 
a 10% mass gap there is a large portion of the parameter 
space in which the WMAP data cannot be satisfied, while 
the regions showed in the second and third plot in fig. 1 
are similar and thus are mass-gap independent. 

3.3. Wino-higgsino NLSP 

If the NLSP is mostly a wino-higgsino we have a three 
particle coannihilation, because the lightest wino is almost 
degenerate in mass with the lightest chargino, so they both 
contribute to the coannihilations. 

In this case we perform the same numerical calculation 
illustrated in the previous subsection. We have extensively 
studied a sample model with mass gap 10%, showing an 
example of funnel region, a resonance that occurs when 
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Figure 1: Plot of the relic density of the LSP vs the stiickelino and 
primeino masses. The first plot shows the case of mass gap 5 %. 
The second plot is a zoom of the first in the region between 20-f40 
GeV for stiickelino mass and 500-850 GeV for primeino mass. The 
third plot shows the case of mass gap 10 % in the same region of the 
second plot. All the masses are expressed in GeV 

2 Mlsp ^ Mao- In our sample model Mlsp ^ 330 GeV, 
while Mao ~ 630 GeV and this leads to the relic density 
plot showed in figure 2, with the same conventions used in 
figure 1. So we can state that also in the case of anomalous 
LSP we can have a behaviour similar to that of the MSSM, 
given that the LSP coannihilates with a MSSM NLSP. 

4. Conclusion 

We have modified the DarkSUSY package in the rou- 
tines which calculate the cross section of a given supersym- 
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Figure 2: Example of a funnel region for a wino-Iike NLSP, whose 
mass gap with the LSP is 10% 



metric particle (contained in the folder ^/src/an) adding 
all the new interactions introduced by our anomalous ex- 
tension of the MSSM. We have also written new subrou- 
tines to calculate amplitudes that differ from those already 
contained in DarkSUSY. We have also modified the rou- 
tines that generate the supersymmetric model from the 
inputs, adding the parameters necessary to generate the 
MiAUSSM Q and changing the routines that define the 
model (contained in '-^/src/su) accordingly. Finally we 
have written a main program that lets the user choose if 
he wants to perform the relic density calculation in the 
MSSM or in the MiAUSSM. The code of our version of 
the package is available contacting andrea.mammarella 
@roma2. infn.it. 

These modifications have permitted to extensively numer- 
ically explore the parameters space for an anomalous ex- 
tension of the MSSM without restriction on the neutral 
mixing or on the free parameters of the model. We have 
verified that our model does not lead to any divergence or 
instability. 

We have also found sizable regions in which we can sat- 
isfy the WMAP data for mass gaps which go from 5% to 
beyond 20%. We have studied some specific sets of param- 
eters for the 5% and 10% mass gap regions, showing that 
relatively small changes in the mass gap can produce very 
important changes in the area of the regions which satisfy 
the experimental constraints. 

We have also showed the presence of a funnel region, analo- 
gous to that in the MSSM, in some region of the parameter 
space. 

So we can say that a model with an anomalous LSP can 
satisfy all the current experimental constraints, can show 
a phenomenology similar to that expected from a MSSM 
LSP and can be viable to explain the DM abundance with- 
out any arbitrary constraint on its parameters. 
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